Linear sigma model at finite density in the 1/A^ expansion to next-to-leading order 
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We study relativistic Bose-Einstein condensation at finite density and temperature using the 
finear sigma model in the one-particle-irreducible 1/A'^-expansion. We derive the effective potential 
to next-to-leading (NLO) order and show that it can be renormalized in a temperature-independent 
manner. As a particular application, we study the thermodynamics of the pion gas in the chiral limit 
as well as with explicit symmetry breaking. At nonzero temperature we solve the NLO gap equation 
and show that the results describe the chiral-symmetry-restoring second-order phase transition in 
agreement with general universality arguments. However, due to nontrivial regularization issues, we 
are not able to extend the NLO analysis to nonzero chemical potential. 
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I. INTRODUCTION 

The Lagrangian of quantum chromodynamics (QCD) 
with Nf quark flavors has an extended SU(iV/)i x 
SU(7V/)i{ symmetry in the chiral limit, which is not re- 
spected by the true QCD ground state. The chiral sym- 
metry is spontaneously broken down to the vector sub- 
group SU(A^/)y by quantum effects, and according to 
Goldstone's theorem there are N't — 1 massless excita- 
tions associated with this breaking. For Nf = 2, this 
implies the existence of three Goldstone bosons, which 
are identified with the pions. The fact that these par- 
ticles are not exactly massless is taken into account by 
small nonzero quark masses in the QCD Lagrangian. In 
low-energy QCD phenomenology, one takes advantage of 
the isomorphism between the groups SU(2) x SU(2) and 
SO (4) to describe the pions as well as the sigma by scalar 
field theory with this symmetry. At low temperature, the 
0(4) symmetry is broken down to 0(3) either explicitly 
or spontaneously. In the latter case, one expects it to be 
restored at sufficiently high temperature. Lattice simu- 
lations of QCD suggest that the critical temperature is 
around 170 MeV, depending on the masses of the quarks. 

At zero temperature, the 0(iV) linear sigma model was 
studied a long time ago using the one-particle- irreducible 
(IPI) l/7V-expansion, to leading order (LO) by Coleman 
et al. [ij], and to next-to-leading order (NLO) by Root 2]. 
At finite temperature but vanishing chemical potential, 
the LO thermodynamics was studied by Meyer-Ortmanns 
et al. [3| . Later several attempts to calculate the pressure 
to next-to-leading order have been made, but the ap- 
proximations made were uncontrolled and renormaliza- 
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tion was ignored. Due to the nontrivial renormalization 
issues, the NLO pressure was calculated only recently Q. 
Specifically, one problem was that the effective poten- 
tial seemed to have ultraviolet divergences which were 
temperature-dependent, except at its minimum. Renor- 
malization of the effective potential could therefore be 
carried out with temperature-independent counterterms 
only on shell. 

Sigma models have also been investigated in detail at 
finite temperature using other methods: Optimized per- 
turbation theory fsl], two-particle- irreducible (2PI) effec- 
tive action either in the Hartree or large- A^ approxima- 
tion @, 0, [^, M, [l3l- III these calculations, one typically 
encounters a local gap equation for a mass parameter, 
which is straightforward to solve. Beyond leading order, 
the resulting equations are nonlocal and much harder to 
solve. 

Linear sigma models at finite chemical potential that 
describe relativistic Bose-Einstein condensation, have re- 
ceived interest in recent years due to the possibility 
of e.g. pion and kaon condensation in compact stars 
[111 . Il2l |. Pion condensation and the phase diagram 
of two-flavor QCD have been investigated using chiral 
perturbation theory [H 111 E [S [13, [II, lattice 
QCD [ii,[23, ladder QCD ^, chiral quark model [H, 
Nambu-Jona-Lasinio models in the mean-field approxi- 
mation [H, [13, M, M, M, [H, [29|, l30|, and the linear 
sigma model using the Hartree approximation [3l|, or 
the large- Af approximation |33 |. 

_ In the present paper, we extend the calculation of Ref. 
[j] using the IPI-l/A'^ expansion by including nonzero 
chemical potential. In contrast to the previous work we 
show that, when interpreted properly, the effective poten- 
tial can be renormalized in a temperature-independent 
manner even off its minimum. Wc apply the method to 
the study of pion condensation. The paper is organized 
as follows. In Sec. [Ill we introduce the general model 
and briefly discuss interacting bosons at finite temper- 
ature and density. In Sec. IIIIl we review the standard 



calculation of the LO effective potential and determine 
the phase diagram. In Sec. HVl we calculate the NLO ef- 
fective potential and discuss in detail its renormalization. 
In particular, we demonstrate how it is used to derive a 
renormalized gap equation. In Sec.|V]we summarize and 
conclude. 



Another virtue of introducing the auxiliary field is that 
the Lagrangian now depends linearly on the two param- 
eters, /^ J and l/Afc. As will be shown later, the di- 
vergences that appear in the effective potential may be 
removed by their suitable redefinition. For the sake of 
systematic expansion in powers of 1/A^, we write them 
as 



II. EFFECTIVE ACTION IN 1/N EXPANSION 

A. The model 

The Euclidean Lagrangian for a Bose gas of N real 
scalars is given by [j| 



^4(«,«' + ^( 



i4>i-Nny--jNH4,n, (1) 



where i = 1,...,N. The parameters /jr.b and At, de- 
note the bare pion decay constant and coupling, respec- 
tively, li H — Q, the Lagrangian has an O(A^) symmetry, 
while if iJ 7^ 0, this symmetry is explicitly broken down 
to 0(A^ ~ !)■ Each generator of the symmetry group 
gives rise to a conserved charge and for 0{N) there are 
N{N — l)/2 such charges. We can characterize a sys- 
tem described by Eq. ([1]) by the expectation values of 
the conserved charges Qi and for each of them we can in 
principle introduce a chemical potential /z^. Nevertheless, 
it is only possible to simultaneously specify two different 
charges if they commute. For the 0{N) symmetry the 
maximum number of commuting charges is | N/2 \ (here 
[ J denotes the integer part of a real number) [33J . 

If we treat two real scalar fields, c/jr and 0/, as the real 
and imaginary part of a single complex field, $ = ((J)e. + 
icj)i)/\/2, the conserved charge Q can be incorporated by 
making the substitution in the Euclidean Lagrangian, 



9o$ ^ {do - m)«', 
do<^^ ^ {do + fi)<^l 



(2) 



The actual value of the charge density in thermodynamic 
equilibrium is then determined from the free energy den- 
sity JF as Q = —dT/d^jL. 

In order to eliminate the quartic term from the La- 
grangian ID), we introduce an auxiliary field a and add 
to Eq. Ill) the following term [l| , 
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The Lagrangian can now be written as 



c = ^{d,ci^,r~'^a{ 






NHcj) 



N- 



(3) 



By using the equation of motion for a, one can elimi- 
nate this field altogether and thus recover the original 
Lagrangian |T])- 



fTT.b 



L+ao + —ai 
A;=A+^"+iV^^^ 



(4) 



where f-^ and A are the renormalized decay constant and 
coupling. The parameters ao.i and 5o,i are the countcrt- 
erms needed to cancel the divergences at the correspond- 
ing orders in 1/A^. 



B. Application to pion gas 

We next use our model to describe the pion gas at finite 
temperature and chemical potential. The vacuum expec- 
tation value (po of (J)n is interpreted as the usual chiral 
condensate. Using Eq. |2]) we introduce isospin chemical 
potential fij for the pair 0i , 02 , which are hence iden- 
tified with the real and imaginary parts of the charged 
pion field. The neutral pion is 03. The Lagrangian |3]) 
thus becomes 



(5,,0.)' - 



,0,-7V/JJ+2^a2 



NH(f> 



N 



if^i{(t>ido(t>2 - 029o0l) - 7rAi?(0i + 02)- 



The Lagrangian is quadratic in the scalar fields (f>i and 
we next integrate over the iV — 3 fields 03, ... , 4>n-i- 
This gives an effective action for the remaining fields 

a, 01,02,0Af, 



1 



5'cff=7r(A^-3)Trlog(-9^ 



la) 



P 



dr / d-^x 



{d,Ar 



{d,hf 



{d,4>Nf 



NH(j)N 



i^/(0iao02 - 029o0l) - ^mK'/'i + (t^l) 
N 7 



(5) 



Now we introduce a nonzero expectation value po for 0i 
to allow for a nonzero pion condensate. We also expand 
the auxiliary field a around its expectation value iM'^. 
Thus we make the following replacements in Eq. ([5]), 



01 -^ VNpo + 01, 
6jv -^ \/iV0o -I- 0Ar, 

f2 



(6) 



a 



iM^ 



a 
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The various factors of v TV do not change the physics, but 
merely facihtate the calculations in powers of 1/A^. Mak- 

I 



ing the substitutions ([6]) in Eq. (O, the effective action 
becomes 1441 



5'eff=-(A^-3)Trlog(-a^ + Af 



/r2 






dr i (Tx 



\{d,cj^if + \{d,hf + \{d,cj^Nf 



Z Zai, Zai, Z 



la (pocpN + Po0i + 



2VN 



(7) 



Expanding Eq. ([7]) up to second order in 1/vN and using Eq. ([4]), we obtain 



|p.i».,.S..2-/;,-^ 



NHcjyo - \N^Ji]pl + \{N- 3)^ ln(p2 + m^) - 1 



U.^i-p) 



I in(p, M) + J + 6, 



-«0o 




p2 + M^ 

2 I ^2 



\ 






- ^NM^ao - -A^Af^foo 



x(P)-iAf2ai-iA/%, (8) 
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where x(^) = ia{P),(f>N{P),MP),MP)f and m^ = 
M^ — /ij . The self-energy function is 



"<''■"' = ^oO^ 



1 



,Q ^- ^ M2 (P + Q)2 + M2 ■ 
We have also introduced the notation 

w„=2;raT' 



t^r^-J- 



^2tt) 



The sum-integral involves a summation over Matsubara 
frequencies and an integral over three-momenta p. At 
zero temperature, the sum-integrals reduce to a four- 
dimensional integral, where we write 



d^P 



(27r)4 

Such integrals are regularized by a four-dimensional ul- 
traviolet cutoff A. 

In deriving Eq. ([8]) we implicitly assumed that the 
classical fields (t)Q,po,M^ are constant. This makes it 
possible to define the effective potential by dividing by 
the space volume and inverse temperature. Specifically, 
we define its LO and NLO parts by 



S, 



off 



pv 



NVlo + V^i^o- 



In equilibrium, the effective potential is then equal to the 
pressure. 

At this point it is worthwhile to discuss explicitly the 
symmetries of the model with(out) the H term and the 



chemical potential, and their spontaneous breaking by 
the condensates (f>0! Po, in the physical case iV = 4. First, 
in the chiral limit, H = 0, the system has 0(4) symmetry 
in the vacuum which is broken explicitly by the chemi- 
cal potential down to 0(2) x 0(2). The pion condensate 
(which is always favored in the chiral limit, see Sec. IIIIBp 
breaks this down to 0(2). The Goldstone theorem then 
predicts a single exactly massless Goldstone boson, which 
is identified with the charged pion (its antiparticle being 
massive). Only at zero chemical potential is the unbro- 
ken symmetry larger, 0(3), leading to three Goldstone 
bosons. 



Second, nonzero H breaks the symmetry of the action 
explicitly down to 0(3). The chiral condensate is a sin- 
glet of this unbroken symmetry, and therefore does not 
affect the low-energy spectrum. Nonzero isospin chem- 
ical potential breaks the 0(3) explicitly down to 0(2), 
which is subsequently broken spontaneously by the pion 
condensate, so we again find a single massless Goldstone 
boson, to be identified with the charged pion. 



Let us remark that, as is clear from the second line of 
Eq. ([5]) , the above discussion of the excitation spectrum 
applies only after the NLO correction to the effective ac- 
tion, and thus the contribution of the chemical potential 
to the kinetic term, has been included. It can then be 
shown using the NLO gap equation that the expected 
Goldstone bosons are indeed exactly massless |J|. 



III. LEADING ORDER 
A. Effective potential and gap equations 

The LO part of the effective potential is given by the 
first fine of Eq. §}), 



V^o = ^M^if. 



00 - Po) + 
i^ ln(p2 + M^) 



2 



H0O 



2 2 

Ai/Po 



A/^flo + -Af'^foo- 



The sum-integral which appears here is evaluated as fol- 
lows. First we extract its zero-temperature part, i.e., the 
four-dimensional integral, which is regulated using a four- 
dimensional cutoff. The difference of the integral over pq 
and the Matsubara sum is UV-finite and calculated e.g. 
using standard contour techniques. The result is 



i 



ln(p2 + m2 



2a2 



167r2 



^ 2 In — ^ 

647r2 I Af 2 



327r' 



:MPM), (9) 



where we have dropped a (quartically divergent) constant 
piece as well as terms suppressed by inverse powers of the 
cutoff, and defined 



MPM) = 



32 

3T4 



P 



dp — n(ijjp). 



Here ujp — ^Jp^^^PhP- , and n(x) = 1/(6^"^ — 1) is the 
Bose-Einstein distribution function. 

The divergences encountered in this calculation can 
be cancelled by adjusting the LO counterterms ao,6o- 
We employ the "minimal subtraction" strategy, i.e., just 
cancel the divergences in the effective potential, leaving 
the finite part unchanged. This leads to 



t^ 



ao 



167r2' 



6o = 



1 



327r2 



1 ^' 



where /i is the scale at which the coupling is renormal- 
ized. The renormalized LO effective potential can then 
be written as 

Vlo = \m\!1 -<g^ pD + ^ Jo(/3Af ) -f i70o 

M^ (•yiTx'^ ^ I? i\ 1 . , , , 

+ 64^(^+^^i?^+2J + 2^?^"- (^«) 

Once we have renormalized the effective potential, the 
gap equations are obtained by setting its derivatives with 
respect to the variables (po, po, and Af ^ equal to zero. 



d(j>o 


opo 


(11) 


^^ -0 




(12) 



As will be explained in detail in Sec. IIVBI it is Eq. (fTTj) 
that represents the true dynamic gap equations, while 
Eq. (fT2|) can be treated as merely a constraint which 
serves to eliminate the vacuum expectation value of the 
auxiliary field a in favor of the condensates (/)o, po- Using 
now the explicit expression for the LO effective potential 
(fTOl). we find 



f 



i)oM^ = H, po{M^ - pj) = 0, 

y2 



pI 



KP_ /32^ 
167r2 (^ A 



In 



Af2 



(13) 



0, (14) 



where 




Ji(/3Af) = 


T^d.h{pM) 8 r 
2 dAf 2 T2 X 




B. Phase diagram 




1. Chiral limit 



dp — n{ujp 

UJp 



In the chiral limit and at zero chemical potential pi, 
the Lagrangian is exactly 0(4)-invariant and the (Jjq and 
Po condensates connected by a symmetry transformation. 
Since switching on the chemical potential favors the pion 
condensate po, we can always set 0o = in this case. 
The second gap equation in Eq. (fT3| then tells us that 
M'^ = pf. Plugging this into Eq. (fH|) . we realize that it 
yields an explicit expression for the pion condensate in 
the chiral limit. 



Po = /' - 



I6n' 



:Ji{Ppi) 



J± 
167r2 



327r^ 



1 ^'' 
m^T 

Pi 



(15) 



At Pi — 0, we use the fact that Ji(0 
this expression to the simple formula 



|7r^ to reduce 



Po 



f^ 



Ji2 

12' 



which predicts a second-order chiral-symmetry-restoring 
phase transition at Tc = \/l2/7r. At nonzero chemical 
potential, the critical temperature cannot be calculated 
analytically. However, at weak coupling it will be pre- 
sumably much larger than pi so that we can still approxi- 
mate Ji{/3pi) with Itt^. We then obtain a weak-coupling 
result for the critical temperature. 



rp2 



12(/^ + ^ 



which is consistent with calculations using different meth- 
ods |32l . |34| . and justifies the assumption we made, 
Tc > pi- 



2. Physical point 

While in the chiral limit the pion condensate persists 
down to zero temperature and chemical potential, i.e., 
the vacuum, when nonzero H breaks the symmetry ex- 
plicitly, the situation changes. The chiral condensate is 
now always nonzero, (J>q = H/M^ , and Eq. p^ modifies 
to 



400 



pI 



f 



IP 



IGtt' 



:JliPfJ'I 



J± 
167r2 



32^ 



In 



This shows that, as expected, the pion condensate in- 
creases with chemical potential and decreases with tem- 
perature. However, at /z/ ^ the right-hand side of this 
equation goes to — oo. That means, there is a critical 
chemical potential fj.jc, given by the solution of 



f 



El 



rp2 

167r2 



2 
Pic 

167r2 



(16) 

below which there is no pion condensate; this is the nor- 
mal phase. In that case AP is found from the gap equa- 
tion (fT4|) by substituting for ^o from Eq. (fT3|l. 



J TT 



El 



rp2 

IGtt' 



(17) 
Note that the normal-phase solution is completely in- 
dependent of the chemical potential. Remembering that 
at the leading order, the parameter M has the meaning 
of the pion mass in the vacuum [4] , the observation that 
Eqs. (fT6| and \n\ are identical upon the replacement 
/i/ ^-> M immediately leads to the physically expected 
conclusion that the critical chemical potential for pion 
condensation is equal to the pion mass in the vacuum. 
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FIG. 1: Phase diagram for pion condensation at the leading 
order. The upper curve is the chiral limit and the lower curve 
corresponds to the physical point. 
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FIG. 2: Pion condensate in the chiral limit as a function of 
temperature and chemical potential. 



3. Numerical results 

We start with discussion of the choice of the parame- 
ters in the Lagrangian ([1]). First of all, in order for the 
theory to be stable, the bare couphng Af, has to be posi- 
tive. This implies that the maximum value the cutoff A 
may take at fixed renormalized coupling A, is 



Amax = M exp 



levr^ 



In Refs. 0, [33, it was shown that the mass of the sigma 
predicted by the model is maximal at A(/i = 100 MeV) ~ 
80, in which case mo- = 433 MeV, lying rather low in 
the experimentally observed range rua — 400 — 800 MeV. 
However, for this large coupling the cutoff is forced to 
be unphysically small since Amax — 720 MeV. We there- 
fore choose a lower value, A(/i — 100 MeV) — 30. The 
symmetry-breaking parameter H is adjusted to H ~ 
(104 MeV)'^ so that it, together with the pion decay 
constant f^^ = 47 MeV (note that our definition of f-^ 



differs by a factor of 1/2 from the standard value), re- 
produces the average of the measured masses of 7r'^,7r , 
m^ = 138 MeV. 

Note that for these values of the parameters, the sigma 
mass comes out as rria — 350 MeV. The discrepancy 
with the measured value might be due to the fact that 
we restrict our treatment to the pion sector, and therefore 
miss contributions from the third- flavor physics. On the 
other hand, it might as well be just a shortcoming of the 
present approach. 

With the above remark in mind, we plot in Fig. [T]the 
LO phase diagram both in the chiral limit and for the 
physical value of H. This phase diagram is identical to 
that found previously in Ref. [33 using the 2PI effec- 
tive action, because the leading orders of the IPI and 
2PI formalisms coincide. In Figs. [2] and [3] we show the 
dependence of the pion condensate on temperature and 
chemical potential, in the chiral limit and at the physical 
point, respectively. 
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FIG. 3: Pion condensate at the physical point as a function 
of temperature and chemical potential. 



IV. NEXT-TO-LEADING ORDER 

A. Effective potential 

The NLO contribution to the effective potential is 
given by the second Hne of Eq. ([8]). By performing the 
Gaussian integral over the field %, we obtain 

^NLO = -^rf ln[/(P,Af)] + ^M^i + ^M%i 

(18) 

Here I{P, M) is the inverse propagator of the a-field in 
momentum space and reads 



/(P, M) = in(P, M) + i + 5o + p42_ 

pliP^+Tn^) 



(P2 _^ „^2)2 ^_ 4^2^; 



2, ,,2 



(19) 



The last two terms of Eq. (|T8|) represent the correction 
to the pressure of a free gas of charged scalar bosons, 
induced by finite chemical potential. It stems from the 
fact that at the leading order, ^i and (/)2 were treated on 



the same footing as all other fields, despite the fact that 
they are endowed with the chemical potential. The finite- 
temperature part of this correction (i.e., the difference 
of the Matsubara sum and the corresponding frequency 
integral) reduces to 



rp4 

647r2 



[K+{PM, P^ii) + K^ipM, p^ii) - 2Jo(/3M)] 



where the functions Kq{PAI, /3fj.i) are the finite- /x/ gen- 
eralizations of Jo(/3M), 

On the other hand, the T = part of the last two terms of 
Eq. (jlSp gives, when evaluated with a four-dimensional 
cutoff, the contribution 



1 



327r2 



^,jK'-^,jM^ + ^ 



This is clearly unphysical since it gives, among others, a 
/i/-dependent quadratic divergence. Should we use un- 
constrained integration over the frequency in combina- 
tion with three-dimensional cutoff or dimensional regu- 
larization in 3 — e space dimensions, we would get simply 
zero instead of these unpleasant terms [36|. However, 
with these regularization schemes, which do not bound 
the size of the four-momentum, we would run into trou- 
ble with the self-energy n(P, M): As we will see shortly, 
it would turn negative at large enough P^ , due to the ex- 
istence of the Landau pole in the running coupling. This 
clash of low-energy and high-energy physics prevents us 
from applying our method directly to the thermodynam- 
ics at nonzero chemical potential [4^. Therefore, from 
now on we set /x/ = and investigate the NLO correc- 
tions to the thermodynamics of chiral symmetry break- 
ing. Moreover, at zero chemical potential we can, without 
lack of generality, set the pion condensate to zero, po = 0- 
The expression for the NLO effective potential (fTS]) 
thus reduces just to the first term plus the counterterms. 
We first evaluate the self-energy function H{P, M), which 
is involved in the effective inverse propagator of the a 
field. Summing over Matsubara frequencies and integrat- 
ing over three-momentum, we obtain 



J 



n(p, M) 



1 



1 A2 ^ 

^ , in — 7T + 1 



P2 + 4M2 , VP2 + 4M2 - VP2 \ 

—^ In 1^====^^= ] + Ut{P, M), 



VP2 + 4Af 2 + ^/p2 ^ 
where IIt{P, M) is the temperature-dependent part of the self-energy, 

1 /•°°, 9, f q^+pq + A^ \ 2 p4 + 4A/2p2 



nT(P,Af) 



87r2p 



4P2 



r 



(20) 



The divergence in the expression P^ coming from the self-energy (PD|) is cancelled by the LO counterterm 6o- 



Proceeding with the same steps as in Ref. |j| , we next 
expand the self-energy in powers of 1/P^ in order to ex- 
tract the dominant ultraviolet (UV) contributions. Af- 
ter averaging over angles, the UV-divergent part of the 
NLO effective potential stems from the part of the ol- 
propagator that we call /uv (-P, Af ) , 



ln[/uv(P,Af)]=lnCi 



where 



1 C. 



1 cl 



P2 Cx 2F4 d 






C2 =2[G - 2Af2) + 2Af^Ci, 

C3 = - 2Af2(G' - ^) - 2M4Ci, 

and we have defined 

K=l + 327r2 ( V + &0 I , 



A 

2 / 2 I rj^l. 



G = 16^"0^ + T'JiipM) - A/2 In ^ ^ . 



In order to avoid infrared (IR) singularity in the 1/ P"^ 
term, we make the replacement P'' — > (P^ + S'^Y ^'^ the 
denominator, where (5 is a small IR cutoff. Needless to say 
that this is just an artifact of the l/P^ expansion and the 
full expression for the NLO effective potential, of course, 
does not depend on 5. We next integrate ln[/uv(^, Af)] 
over P and extract the UV-divergent part. The result 
can be written as 

- i /" ln[/uv(-P, M)] =D + UV-finite pieces, (21) 
2 Jp 

where the divergent terms are contained in 



D 



p2 



M^ 



G - 2AP 



2M* 



In^i; (P2+J2)2 



167r2 



-A/^A^ + (G - 2Af ^jA^e'^ h ■ 



-2Af4 In 



A: 

5^ 



(22) 



We have provided an exact representation of the diver- 
gence D, which is a necessity for the numerical imple- 
mentation of renormalization. 



B. Renormalization 



The UV-divergent term in Eq. ([H]), proportional to G, 
depends explicitly on temperature, except at the solution 
of the LO gap equation ((14]), where G reduces to IGir^f^. 
This led previously to the conclusion that the NLO 
effective potential is renormalizable in a temperature- 
independent manner only at the minimum [j|. At this 
point, we should carefully distinguish what effective po- 
tential we are speaking of. It is crucial to note that what 
we have actually calculated so far, is not an effective po- 
tential of just the scalar fields 4>i, but one of (j)i and a. 

As is clear already from the LO expression (|10p . such 
an effective potential does not have the usual physical 
interpretation as the minimum energy at fixed vacuum 
expectation values of the fields of the theory. In fact, 
it is unbounded from both below and above and the so- 
lution of the LO gap equations P^ . dH]) is its saddle 
point rather than a global extremum. This problem can 
be traced back to the fact that the field a is not an in- 
dependent dynamical degree of freedom; upon using its 
equation of motion, it turns out to be a composite oper- 
ator function of (pi . 

The way out is to treat the gap equation for M p^ as 
merely a constraint that serves to eliminate M in favor of 
(f>Q. We thus arrive at an effective potential as a function 
solely of (f>Q, which already has the usual properties [l|. 
This is the physical effective potential, V, we eventually 
want to calculate, independent of the method we use. 

First of all, let us note that none of the results of 
Sec. mil is altered by changing the perspective from 
VLo(Af, (/)o) to Vlo(0o)- The reason is that, of course, 
the stationary point of Vlo coincides with that of Vlo 
once the LO constraint (|14p has been used to eliminate 
M from it. 

At the next-to-leading order we have to be more care- 
ful. Expanding the solution of the constraint P^ in pow- 
ers of 1/A'^ (and suppressing the argument in A/((/)o) for 
the sake of legibility) , the physical effective potential up 
to next-to-leading order reads, 



V(0o) = iVVio(A/Lo + ^A/nlo + • • • , 0o) + VNLo(AfLO + ^A/nlo + ■■■ ,M 



NVw{Mw, M + ?^(AfLO, 0o)AfNLO + Vnlo(A/lo, 0o) + 0{1/N) 
oM 



NVlo{Mlo, 0o) + ■t^NLo(AfLO, M + 0{1/N), 



where the middle term vanishes due to the LO constraint (fT4ll . We can thus see that in order to calculate the 
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physical effective potential to next-to-leading order, it is 
sufficient to solve the LO constraint for M. Doing that, 
we obtain the effective potential as a function of (j)o alone. 
Moreover, Eq. (P^ tells us that the divergences of Vlo 
are independent of temperature for all values of (Jjq since 
G = IGn^f^ follows directly from the constraint (fT4|) and 
does not need the true gap equation (fT3|) to be fulfilled. 
Substituting for G, Eqs. ^ and HH]) show that the 
necessary NLO counterterms are 



A2 
ai = -^ ( 1 -K 2e" h ^ 



1 A2 
6i = --^ln— . (23) 



jj, the effective potential 



Finally, it is important to check whether the procedure 
of elimination of M is in practice well defined. To this 
end, note that the left-hand side of Eq. (fT4|) is a mono- 
tonically increasing function of M^. It goes to 4-oo as 
IVP — > -f oo, and reaches the limit f^—cpQ — j^ as M^ -^ 0. 
This means that the constraint has a real solution only 
when (/)Q is such that (J)q > f^ — j^- If this condition 
is satisfied, the solution of the constraint is unique as it 
should. Inside the interval {—R, +R) around the origin 

with the "radius" R = JJl 

V(0o) is undetermined by the auxiliary field technique. 
In principle it can be analytically continued to this inter- 
val [371 (this is what we actually do in our weak-coupling 
analysis in Sec. IIVC|) . but that would be hard to imple- 
ment numerically. Fortunately, for our parameter set the 
minimum of the NLO effective potential always lies out- 
side this interval so that the problem does not arise. It 
is also useful to note that the interval (— i?, +R) is pre- 
cisely the region where the scalar field fiuctuations be- 
come tachyonic, and the perturbative one-loop effective 
potential acquires nonzero imaginary part. 



C. Gap equations 

Once we have renormalized the NLO effective poten- 
tial, we may proceed in the standard manner and calcu- 
late the condensate (p^ as a solution to the gap equation. 



dV 



0. 



(24) 



Thanks to the fact that M now depends on the conden- 
sate, we have to evaluate the total derivative using the 
chain rule, which leads to 



dVi^o l_ / 9Vnlo , S^NLO dMl 



LO 



dAP 



0. 



(25) 



It is important to observe that while both 9Vnlo/90o 
and dVjsiho/dM^ may in general contain temperature- 
dependent divergences (just because we can set G = 
IGiT^f^ only after taking the derivative), these cancel in 
the linear combination in Eq. (j25p . 

To illustrate this procedure by an explicit example, 
let us consider the weak-coupling limit in case of exact 



chiral symmetry, i.e., H = 0. Since Af = in the chirally 
broken phase at the stationary point of the LO effective 
potential (fTO|) . we expect it to be of order 0{X) at the 
NLO value of the chiral condensate. Then, in the LO 
constraint p4)) we can set M — everywhere except the 
term proportional to 1/A, which leads to the expression 
for M, 



M^ = 



A 



-f 



12 



(26) 



The first partial derivative in Eq. ([25]) is simply —M(po, 
while the last one is A^o by Eq. (|26|l. The 0{\) part of 
the derivative of Vnlo with respect to (J)q is given by 



dVf 



NLO 



d4'Q 



= -(Po 



1 



p2 + M^ 



in(P,Af) + i + 6o + 74^ 



-A0, 



4pi^ 



M2- 



By an analogous argument, the partial derivative with 
respect to AP is ^ai + 0(A) -^ A'^/lGir^ as A ^ 0. 
Putting all the pieces together, we find 



dV 



-M^o 



N' 



'° (^P P^ 



A2 



AP 



AP 167r2 
N 12 



O(A'). 



The sum-integral above was calculated by differentiation 
of Eq. ^ with respect to AP. This, together with the 
LO constraint ((26| . immediately yields the NLO expres- 
sion for the chiral condensate in the weak-coupling limit. 



J TT 



2\ T^ 



and the NLO critical temperature 



T, = 



12 



1 



h- 



(27) 



It is important to stress that while the LO critical tem- 
perature, Tc = vl2/,r, derived in Sec. IIIIBl is a nonper- 
turbative result valid for all values of the coupling, the 
present NLO correction is just a weak-coupling limit. In 
fact, while Eq. (j27[) predicts a decrease of the critical 
temperature with respect to the LO value, in particular 
for iV = 4 by about 20%, our numerical analysis dis- 
cussed in Sec. HVP] shows that at A(^ = 100 MeV) = 30 
the critical temperature actually increases by more than 
30%! In order to make sure that these two results are 
consistent, we will briefly investigate the dependence of 
the critical temperature on the coupling. 

Apart from calculating the exact minimum of the NLO 
effective potential, the NLO correction to the chiral con- 
densate can also be determined in a more straightfor- 
ward manner [38[. One writes the solution to Eq. 



as (j)o = 0Q° + j^4'o^'^ ^^^ expands the whole gap equa- 
tion in powers of 1/A^. Comparing the terms of the same 
order, one thus gets an exphcit formula for the NLO cor- 
rection to the condensate, 



60 



iNLO 



rfVNLo(0o°) 

d<t>o 

d^VLoi4>};") 

d<Pl 



(28) 



While the first derivative of the NLO effective potential in 
the numerator of Eq. ([28)1 is best evaluated numerically, 
the second derivative of the LO effective potential in the 
denominator can be calculated analytically. Keeping in 
mind that this is a total derivative we find, in a fashion 
similar to Eq. ([25]) . 



#0 9(j)l 



d^Vj 



LO 



9(Af2) 



2^2 






At the intermediate step, we used the implicit function 
theorem to extract the derivative dMlQ/d4>o from the 
constraint (fT2|) . The second partial derivatives are easily 
found to be d^VLo/d(l)l = ~A'P, d^VLo/d(j)odAP = -(j)o, 
and 



ao 



1 



a(M2)2 

where 



327r2 



J2(/3M) 



32^2 



In 



M2 



1 



(29) 



J2(/3M) = -T 



dNP 



= 4 



dp 



n{ujp) 



OJv 



D. Numerical results 

For sake of numerical computation, the NLO effec- 
tive potential was renormalized as follows. First, the 
"vacuum" effective potential, evaluated at zero temper- 
ature and with the LO value of the chiral condensate, 
was subtracted. This canceled the leading, quartic, field- 
independent divergence. In addition, we subtracted the 
integral representation of the remaining divergence (j22p . 
This amounts to NLO rcnormalization according to Eq. 
(P5)) as well as a finite piece, which must be added to the 
effective potential after the numerical integration. 






M4 



In 



,52 



where Mq refers to the value of the mass parameter in 
the vacuum. 

Note that the UV-finite pieces indicated in Eq. (PT|) 
depend explicitly on the counterterm bg, hence also on 
the cutoff. It is technically advantageous to keep these 
terms and regulate the integrations with the cutoff even 
after all divergences have been subtracted. Physically 
this means that we treat our model as a low-energy ef- 
fective theory and neglect contributions suppressed by 
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FIG. 4: Solution of the NLO gap equation in the chiral limit. 
Upper panel: the chiral condensate; lower panel: the mass 
parameter M. For reference, the LO values are shown by the 
dashed lines. 



inverse powers of the cutoff. In Ref. [J] the cutoff depen- 
dence of the final results was shown to be mild enough to 
justify this procedure. In concrete calculations, we chose 
A = 5 GeV. 

In Fig. |4]we display the numerical solution of the NLO 
gap equation (|25p in the chiral limit. We can see that the 
correction to the chiral condensate is quite large. How- 
ever, at zero temperature, the increase is about 25%, 
which is exactly what one would expect from a 1/iV cor- 
rection at A^ = 4. The critical temperature increases 
from the LO value 163 MeV to 218 MeV, which means 
an increase by about a third. That is still compatible 
with the 1/A^ prediction up to a numerical factor close 
to one. 

At the physical point, the NLO corrections are con- 
siderably smaller, see Fig. \E\ Here we also compare the 
exact self-consistent solution to the NLO gap equation 
(solid line) with the expansion of the condensate accord- 
ing to Eq. (PSJ) (dash-dotted fine). [In the chiral hmit 
it cannot be used since the function J2(J3M) diverges 
at the LO solution A/lo = 0.] Apparently, the formula 
(PS)) . though exact in the large- A^ limit, is of little use in 
our case: At temperatures below about 200 MeV, it over- 
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FIG. 5: Solution of the NLO gap equation at the physical 
point (sohd hues). Upper paneh the chiral condensate; lower 
panel: the mass parameter M. For reference, the LO values 
are shown by the dashed lines. We also show the NLO chiral 
condensate calculated using Eqs. (|28l) (dash-dotted) and (|30p 
(dotted), respectively. 



shoots the NLO correction to the condensate provided by 
exact solution to the gap equation (P5)l by nearly 100%. 
In spite of that, the two evaluations of the NLO conden- 
sate are formally consistent for their difference is of order 
l/N"^ of the LO value. 

However, the agreement with the exact solution of the 
NLO gap equation can be improved by modifying Eq. 
(1281) to 



rfVNLo(0g^ 



iNLO 



Xa2 



N 






(30) 



This expression is easily understood. The denominator 
contains the second derivative of the full (LO plus NLO) 
effective potential, and the numerator contains the first 
derivative of the same, just because dVLo(0o'^)/'^0o = 
0. Eq. ((50)) then gives the approximate solution to the 
NLO gap equation by replacing the function dV/d(j>o at 
the point (J>q = ^q*-* with a straight line; this is just the 
first step of the Newton iterative algorithm for solving 
nonlinear equations. 



200 
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FIG. 6: Critical temperature for chiral symmetry breaking 
as a function of the coupling. The dashed line denotes the 
LO critical temperature while the dotted one indicates the 
weak-coupling NLO limit (|27p . 



From the point of view of 1/N expansion, Eq. ([301) 
resums a selected series of contributions to the gap from 
arbitrarily high orders. Our comparison to exact solu- 
tion of the NLO gap equation demonstrates that Eq. 
(|30)) (dotted line in Fig. [5]) is superior to the simpler 
expression (|28p at a little extra computational cost (just 
one more evaluation of the NLO effective potential when 
implemented properly). 

In order to make connection between our numerical re- 
sults and the weak-coupling prediction (j27p for the crit- 
ical temperature in the chiral limit, we also investigated 
the dependence of the critical temperature on the renor- 
malized coupling, with all other parameters fixed. The 
result is shown in Fig. [HI We could not check the limit 
(PT)) directly, because as follows from the discussion in 
Sec. IIVBI the auxiliary field technique only allows us to 
study the range of coupling for which the NLO critical 
temperature is higher than the LO one (and the chiral 
condensate at fixed temperature thus larger). Neverthe- 
less, the plot in Fig. [H] clearly indicates that the critical 
temperature decreases at small coupling and the corre- 
sponding numerical values are consistent with the limit 
(P7)) . (The bend at large coupling is most likely a cutoff 
effect: At A = 40 the maximal cutoff allowed by stability 
criterion is already very close to the actual value used in 
the computations.) 

In addition to the nonuniversal quantities such as the 
critical temperature and the condensate 0o a-s a func- 
tion of temperature, our NLO calculation also allows to 
check some universal properties of the 0{N) model. At 
high temperature, the nonzero Matsubara modes decou- 
ple and the system undergoes dimensional reduction to 
a Euclidean field theory in three dimensions. This im- 
plies that the critical exponents of the system are those 
of the three-dimensional O(A^) model. For example, the 
critical exponent i^ that governs the behavior of the or- 
der parameter 00 near the critical temperature, has been 
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calculated in the 1/A^ expansion to next-to-leading order 

MM, 



1 



1^ = ^ — 



Ntt'^ 



©(l/iV^). 



For A^ = 4, this reduces to z^ w 0.3987 (as compared to 
the LO value v = 1/2). The power-law fit based on the 
NLO chiral condensate in Fig. 3] is in agreement with 
this analytical result within the error bars. 



V. CONCLUSIONS 

In the present paper, we have investigated the thermo- 
dynamics of the 0{N) linear sigma model at the next-to- 
leading order of the IPl-l/N expansion. We explained 
how the NLO effective potential can be systematically 
renormalized in a temperature-independent manner. The 
crucial observation leading to this result was that the ef- 
fective action previously used in literature, is one of the 
scalar field 0i as well as the auxiliary composite field a. 
The latter has to be eliminated before renormalization 
can be carried out consistently. One thus arrives at an 
effective action which is a function solely of (pi and has the 
usual physical properties and interpretation. This proce- 
dure is further justified by the fact that in the numerical 
computation, we solved the NLO gap equation by direct 
extremization of the effective potential, thereby proving 
that the found solution provides a thermodynamically 
stable configuration. 

A different way to calculate the NLO corrections to 
the condensates, based on expressions like (^5)) . is some- 
times used in the literature. This has the great advantage 



of being much less computationally demanding than the 
self-consistent solution of the NLO gap equations, in ad- 
dition to the fact that the latter may lead to unphysical 
results [41[. In our model the NLO gap equation can be 
solved exactly. A direct comparison of the results ob- 
tained using the two methods shows that a careless use 
of Eq. (|28p can be misleading and need not improve the 
approximation provided by the leading order. We sug- 
gest a simple modification, Eq. ((30)) . which in our case 
yields much better results with a little extra computa- 
tional effort. 

We also extended the NLO IPI formalism to finite 
chemical potential. At leading order, we determined the 
phase diagram and confirmed the results previously ob- 
tained using the 2PI effective action at the leading order 
[33 | . This is of course not surprising since at the leading 
order, the IPI and 2PI formalisms coincide. It would be 
interesting to compare these two approaches at the next- 
to-leading order. In the present framework, we could not 
extend the explicit NLO calculations to finite chemical 
potential due to regularization problems. This issue as 
well as the relation to the 2P1 formalism [i^, |43| will be 
subject of future work. 
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We are omitting terms linear in the quantum fields as 
these terms vanish at the minimum of the effective po- 
tential. 

Nevertheless, it should be stressed that the strategy for 
extracting the UV-divergent part of the effective poten- 
tial would formally not change even at finite chemical 
potential. The only necessary modification of the follow- 
ing formulas for the NLO effective potential would be 
the natural replacement 0o ^ 0o + Po in Ihs definition 
of the quantity G. The divergent part of NLO effective 
potential (|22p . as well as the whole NLO renormalization 
procedure including the prescription (|23p . would remain 
unchanged. 



